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Abstract 

In the context of the Grand Unified MSSM, we investigate the fermion mass matrices 
at GUT scale. We note that from the experimental mass pattern the determinants of the 
Yukawa matrices at this scale can be unified with good precision. Taking the unification 
o determinants as an hypothesis, it gives two model independent predictions that in the 
MSSM turns out to determine an appropriate value for the product m,im s and tan/3 ~ 7-10 
in the favored range. We then review a predictive model of SU(3) flavour in the context of 
supersymmctric SO(10) that nicely implements this mechanism, while explaining all fermion 
masses and mixings at la level, including neutrino data. 



1 Introduction 

Understanding the physics that underlies the measured pattern of fermion masses and mixings 
is still an open problem. While there is little clue on the way to follow, it is certain that 
one will need a framework beyond the Standard Model. We approach here the problem in 
the interesting setup of SO(10) Supersymmetric Grand Unification (SUSY GUT) with an 
horizontal SU(3) flavour symmetry [2]. The reason for this choice is of course that all the 
fermions of a Standard Model family plus the right-handed neutrino fit in a basic multiplet 
of SO(10), the 16 representation, and that the SU(3) flavour symmetry is the largest allowed 
(actually U(3)) by gauge interactions. SUSY provides a one-scale unification of the gauge 
couplings, in addition to benefits like the stability of the (large) hierarchy, and others. 

While in the low energy world both the gauge SO(10) and flavour SU(3) symmetries are 
broken, in this framework the Yukawa matrices are treated in an unified way, so that they will 
show the remnants of both these breakings. 

With this in mind we approach the problem by studying the mass matrices at the Grand- 
Unification scale Mqut — 2-10 16 GeV. We consider first the charged-fermions Yukawa couplings 
renormalized at this common scale assuming MSSM running (recalled in appendix EJ. We 
plot them in figure Q showing also the ranges resulting from experimental uncertainties, and 
illustrating the three relevant cases of small, moderate and large tan (3. 

The first observation that one can make in figure ^ is that the Yukawa couplings show clear 
hierarchic patterns: indeed if it were not for some deviations, the picture would be quite simple: 
If ~ e| : ef : 1, with small parameters ef, for / = u,d,e. One can estimate the hierarchy 
parameters as c u ~ 1/500 in the "up" sector, while in the "down" and "charged leptons" 
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Figure 1: Yukawa couplings at GUT scale calculated from experimental charged fermions 
masses, assuming MSSM running and including uncertainties also from l± renormalization. U, 
D, E are respectively red (continuous), green (dashed) and blue (dot-dashed). 



sectors the hierarchies are very similar, e d ^ e — 1/50, reaching the well known observation that 
e u < e e ~ e d . 

About the deviations, it is evident that the "up" sector (U) has an almost exact hierarchy, 
while the "down" and "charged leptons" sectors (D, E) deviate sensibly. These deviations can 
be thought to belong mainly to the first two generations, so that the picture that emerges is 
the following: 

l u : l c : l t ~ el : e u : 1 (1) 

Id ■ h ■ h ~ -j-ed : A d e d : 1 (2) 

le ■ lp ■ It ~ ^7-e 2 e : A e e e : 1. (3) 
A e 

The factors A e , A d that parametrize the deviations can also be estimated as A e ~ 2.5 and 
A d ~ 0.5. 

Other noticeable facts that can be deduced qualitatively from figure^are: 1) b-r unification, 
lb — l T , that works in all tan/3 regime; 2) at large tan/3 ~ 50 also "top" is unified, It — lb — It 5 
3) at very small tan (3 ~ 2 the first family is approximately unified, as is the second family at 
intermediate tan/3 ~ 7. Aspects 1) and 2) were analyzed deeply, see e.g. pj], while very small 
tan j3 was suggested as a possible inverted hierarchy unification [UIS1- Here we will focus on the 
unification of the second family, that works at intermediate tan (3 and leads to an interesting 
scenario. 

As far as neutrinos are concerned, we have no direct information on their Yukawa couplings: 
even though we now know that they have a mass and we know their hierarchy and mixing angles, 
the right-handed neutrinos, that are automatically present in SO(10), may have a separate 
Majorana mass matrix. Therefore one should not expect neutrinos to follow the same scheme 
of the charged fermions, and will be considered later. 
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2 Unification of determinants 



A step beyond these observations can be made by using together two known "vertical" relations 
between the D and E sectors, namely that lj, s — lelu, and If, ~ It. 1 Together, these two relations 
lead to the conclusion that for the D, E sectors one has 

le In It ^ Id h k • (4) 

In words, since what appear here are the eigenvalues of the Yukawa matrices, this simple relation 
suggest that in the D and E sectors the determinants of the Yukawa matrices are unified at 
GUT scale. 

It is natural to try to extend the previous argument also to the U sector. In doing so, in the 
MSSM, we encounter the presence of the tan (3 parameter, that measures the U-D Higgs VEV 
orientation and shifts the U sector Yukawa eigenvalues with respect to the D, E ones; this can 
be seen in the three cases of figure ^ 

Since tan (5 is very loosely bounded, until SUSY is discovered, we may proceed in the inverse 
direction: we may assume the unification of the Yukawa determinants, as if it was motivated 
by some symmetry principle, and look for its phenomenological consequences. 

We find the first result in terms of tan/?: if we impose 

le l/i It = Id Is lb = lu le h j (5) 

it follows that among the three regimes of small, moderate or large tan /3, only the middle one is 
selected. Therefore unification of determinants gives a prediction of tan/3 ~ 7. This qualitative 
conclusion may be checked in figure ^ but we will describe below a precise test of it. It is 
interesting to note that this value of tan (5 lies in the currently favored range |0j • 

Then, still reasoning in the inverse direction, after choosing the appropriate tan (3 ~ 7 
that guarantees the unification of the U with the D, E determinants, we have that all Yukawa 
eigenvalues are built around the common Yukawa scale I defined precisely by I s = l e lnl T , since 
the "lepton" Yukawa couplings are known with best accuracy. Now we note that this scale acts 
as a "pivot" for all the hierarchies: since fermions are hierarchic, the first family fermions get 
lighter when the third family ones are heavier. The middle scale can be approximately identified 
with l c , because the U Yukawa couplings are hierarchic to a very good approximation. 

In the U sector we find the first and simpler example: the more It is large the lower lies l u : 

I 2 

lu = T ' with ^ - l el^T • (6) 

h 

One can think of this mechanism as a kind of seesaw between different flavours, i.e. a flavour 
seesaw. 2 

In the D, E sectors similarly, since l T , If, are large then l e l^ and l^h will be small. Moreover, 
a similar phenomenon takes place: using bottom-tau unification, If, ~ l T , one has that Z e L ~ 
Idh = l 3 /lb,T- This relation says that in the D and E sectors, thanks to b-r unification, in 
the first two generations there is a smaller seesaw: the deviations of L and l e are in opposite 
directions with respect to exact hierarchy, and the same happens for the deviations of Id and l s . 

Of course this mechanisms can not explain why It or If, are higher than the common I, like 
the celebrated seesaw mechanism in the neutrino sector can not explain the high scale of the RH 
neutrino Majorana mass. The unification of determinants only provides the two relations © 

The first holds precisely but with 40% experimental uncertainty from mdm s , the second is exact to 10%, 
with m;, known to 5% accuracy. 

2 Of course until a realization in flavour space is given, the flavour-seesaw mechanism described here has not 
the commonly assumed meaning of 'mixing with heavy decoupled states', but it has the original meaning of the 
seesaw: a quantity going down whenever an other goes up. . . 
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that should be accompanied by a complete "theory of fermion masses" to generate the actual 
hierarchies between yukawas, and also their deviations. 

We will describe below one such a model of Flavour SUSY GUT that includes a natural 
realization of these ideas, but before we will analyze the model independent predictions that 
strictly follow from the unification of Yukawa determinants ® alone. 

The conditions (jSJ give two testable relations between quarks and lepton masses. In partic- 
ular, running down to low energy (see appendix EJ: 

m d m s — — ^ B t {() 

m e mfj,m r Rlr] c 7] uds 3 12 

m u —r-^ [tan p ) B t (8 J 

m c m t Rt-qj^rir 

and these translate into the following predictions for m d m s and m u : 

m d m s = 880 MeV 2 ± 2% mt ± 4% mb ± 25% Q3 (9) 
/tan 3 \ 3 

m u = ( — — J 1.95MeV±20% mt ±10% mc ±15% a3 . (10) 

Using a mean value of 19.5 for the known ratio m s /m d , the first equation gives the central 
values 

m d = 6.7MeV, m s = 131MeV, (11) 

that fit well the data, with m s a bit high, within la of the experimental range. This means that 
from experimental data the D, E determinants are unified at GUT scale within 1-a. Lowering 
a^\z helps this unification, so does raising the SUSY scale. 

The prediction for m u carries a dependence on tan /?, that can be used to derive a prediction 
for it, as a function of m u itself: 

tan/? = 7.3 ± 1.3 m „ ± 0.5 mt ± 0.2 mc , (12) 

where we used the la range for m u = (3 db 1.5) MeV [7j. This prediction lies exactly in the 
favored range determined from the recent analysis in various MSSM scenarios jTHH]. 

An alternative prediction for tan (3 may be derived by combining the two above and assuming 
for m u /m d and m s their experimental ranges. One gets: 

tan/3 = ( m " m c m t \ l,Z R n ( VbVuds \ 1/3 B -n/3 
\m d m s m b ) R d \ i] c J * 

= 10 ± 1.3 mu/md ± l ms ± l mt ± 0.5 a3 . (13) 

This prediction tells what is the value of tan/3 once a model has been solved and m u /m d , m s 
and nit are determined. 

We can note however that the dependence on a^\z cancels almost perfectly in both the 
ratios involving the i?'s and the r/'s, and enters only through the Bt dependence (<4%). This 
is a remarkable fact, considering that for individual masses the uncertainty in a s \z generates 
the dominant errors in their predictions. 

In addition, and equally remarkably, the even larger dependence on the supersymmetry 
breaking scale also cancels almost exactly in the ratio R u /R d , and enters just through the Bt 
factor, where it can be estimated as a small decrease of 3% when Msusy is raised from 167 GeV 
to 1 TeV. 
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Figure 2: Renormalization of Yukawa determinants, from Mz to Mqut in the MSSM, assuming 
tan j3 ~ 8. The left plot shows the three Yukawa determinants for U (red, continuous), D (green, 
dashed), E (blue, dot-dashed). The right plot is more quantitative and shows the D/E ratio, 
with the one-sigma uncertainty. It becomes compatible with 1 beyond 10 8 GeV. 



Before closing this section, it is also interesting to discuss and display the RG running of 
the determinants, from Mz to the GUT scale. We remind that the determinants of the Yukawa 
matrices are not RG invariant: for example for tan/3 < 25 they follow the equations: 

4- In det If = G\ on + C f a t , (14) 
at 1 

where t = 4tt ln(/i/GeV), 0^=1,2,3 are the gauge couplings and at = £?/47T. The beta func- 
tion coefficients C u ^ e = {12,1,0} are due to "top" loops and G u = {—13/5,-9,-16}, = 
{-7/5, -9, -16}, G e = {-27/5, -9,0} are due to gauge loops. 3 

One can check analytically that these coefficients allow the unification of the D, E deter- 
minants at high scale, within their uncertainties. For illustration we plot in figure |2] how the 
determinants are renormalized from Mz to the GUT scale, for tan (3 ~ 8.5. In the right plot we 
show the ratio of the D and E determinants, as a test of the unification, given that the U sector 
is used to define tan j3 as in (|12|) . (|13[). We show also the one sigma uncertainties, defined when 
all "down" masses are deviated simultaneously at their one-sigma limit. The picture shows 
indeed an approximate unification of determinants at high energy, since the curve is compatible 
with the value of 1 for scales above 10 8 , to 10 16 GeV, with a slight preference for a scale of 
~ 10 13 GeV. At the GUT scale the unification is precise within 25%, that is less than the 
uncertainty on the D determinant (see indeed ©)• 

We also comment on the results of a similar analysis in the case of the non-supersymmetric 
Standard Model. We find an even better agreement of the D, E determinants (within 2%) at 
the scale of the a\-a% unification, Mqtjt-12 — 10 13 GeV. However, even ignoring the lack of 
simple unification of the three gauge coupling constants, the requirement of the U determinant 
to match with D, E requires an extension to models of the 2HDM type, with the relative tan (3 
parameter again in the moderate range. It would be interesting to check these ideas in such a 
non supersymmetric model that unifies via an intermediate scale. 

3 Model realization: Rank-1 combinations in flavour space 

In this section we describe a model that realizes the above ideas in the framework of SO(10) 
SUSY GUT with SU(3) flavour symmetry. The question to be answered is "what in this frame- 
work generates fermion mass hierarchies and flavour mixing?". Indeed one can use the fact that 

3 For the (non-supersymmetric) Standard Model, C Utd , e = {21/2, 15/2, 9}, while G„ = {-51/20, -27/4, -24}, 
G d = {-3/4, -27/4, -24}, G e = {-27/4, -27/4, 0}. 
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some hierarchic parameters are generated by the SO(10) breaking: they distinguish between 
the U, D, E sectors, but independently of the flavour. Then one can imagine that when flavour 
is broken these hierarchic parameters will be assigned to the correct direction, i.e. large ones in 
the third family direction in all fermion sectors. Since the unification of determinants observed 
in the previous section seems to suggest that some relations among eigenvalues may not depend 
on the fermion sector, and also do not depend on the flavour mixing, we can then try to think 
that in first approximation it is possible to disentangle the hierarchy of eigenvalues from the 
mixing structure, that in turn can be common to the U, D, E sectors. Once a flavour breaking 
is introduced it will not only generate mixings, but, if it is not just unitary matrices, it will also 
modify slightly the eigenvalues, and thus will generate the deviations A e ^. 

Following ,8, we review here how these ideas can be implemented, writing the Yukawa 
matrices at GUT scale as a hierarchic combination of three fixed rank-one projectors in flavour 
space. This approach was devised in [3J [H] as a description of an inverse hierarchy unification 
at very low tan/3, as in the first case of figure ^ It turns out however that it is quite successful 
in accounting for fermion masses and mixings even at moderate tan f3, in fact it is even more 
successful. 

Let us start with the case of no mixings and ignoring the deviations from strict hierarchy. 
As mentioned in the beginning, it is clear that this case is quite simple, and one would build a 
model that produces the following Yukawa matrices: 



7llO 
'J 




f = u,d,e 



with e ei( i ~ 1/50 and e u ~ 1/500. Notice that the determinant here is independent of the 
fermion sector /, i.e. det If = I s , while the eigenvalues follow different hierarchies in different 
sectors. Unification of determinants follows from choosing {ej, 1, ej 1 } coefficients. 4 This form 
ignores the deviations as well as the mixings, but the reason of this is evident: three chosen 
matrices are mutually orthogonal projectors. 

One is therefore led to consider a modification of the above, using the most general rank- 
one projectors, in place of orthogonal ones. It is necessary that they are rank-one, because the 
eigenvalues should reflect the hierarchy of the €f. 

The Yukawa matrices are thus taken as a hierarchic combination of three generic rank-one 
projectors in flavour space -Pi, 2,3: 

If = l[e f Pi+p2+eJ 1 P 3 \, f = u,d,e,u D (15) 

where we have included also the Yukawa matrix for neutrinos, that will anyway be present 
in SO(10) unification, and we will consider it later. The -P n =i,2,3 are three generic rank-one 
complex matrices and may be parametrized as P n = p n ®p' n = PnP'n, with p n , p' n three couples 
of vectors in flavour space (flavour triplets). The p n and p' n may be realized as the VEV of 
flavour triplets, and may be thought as providing the 'left' and 'right' flavour structures. They 
are thus responsible for the breaking of the SU(3) flavour symmetry, and this breaking does not 
depend on the fermion sector. Once both SO(10) and SU(3) are broken, the ej hierarchies are 
mixed via the P n 's to generate both almost hierarchic eigenvalues as well as mixing angles. 

In this way one has a cooperation of the gauge (vertical) and flavour (horizontal) sectors, 
of the theory, while allowing quite general mass matrices. 

4 We anticipate that since these are built as effective operators, there are two other allowed choices that shift 
the exponents by ±1 power of e/, leading to the other cases of figure^ For example {e^, e/, 1} leads to the third 
case. The first and the third case do not unify the determinants but may be analyzed similarly. They however 
are in high tension with the experimental data: they predict mainly low mt; see [S]. 
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Looking at the determinant of these If matrices, it turns out to factorize as: 



detL 



I Dp Dp/ , 



(16) 



with Dp = pi ■ (p2 AP3) and similarly for D p >. This is again independent of €f, i.e. the deter- 
minants are unified in the U, D, E, v sectors and the eigenvalues follow the predictions of the 
previous sections, independently of the parameters entering in p n and p' n . 

To see in detail how this ansatz works, we choose a flavour basis and we restrict for simplicity 
to the case of symmetric projectors for example assuming that p' n = p n are the VEV of the 
same triplet: 



U = I 



e/A +A +eJ 1 P 3 






( 





c ^ 













H 





1 




K 1 ! 















Vo 


c 


h 


Vo 





d 2 J 





be f 
c + abef 



c + abe f d 2 /e f + c 2 + b 2 e f J 



(17) 



where we have rotated globally all the p n vectors in flavour space and rescaled e/, I to set 
pi = (l,a,6), p2 = (0,1, c), P3 = (0,0, d). Eliminating also the irrelevant phases, we end 
up with the following generic 13 parameters domain: I, a, d £ M + , b,c,e u ,ed,e e £ C. These 
parameters should account for all Dirac couplings of charged fermions. 



It is also worth displaying an other form of this matrix 

ft 



h 









T 



-/ 



Pi 



P2 



P3 




(18) 



where we can see that the flavour breaking triplets p n induce a flavour mixing through T that 
is not unitary, and therefore will induce not only mixing angles, but also deviations in the 
eigenvalues. The deviations will therefore be linked to the flavour mixing. This form is also 
useful to easily prove (fTB]) . 

3.1 Leading order analysis 

The form (|17jl can be taken as an ansatz and can be solved quite completely by computing the 
leading order expressions of eigenvalues and mixing angles: 



/ ~ / \f I 



l + a 2 e d \ 



L ~ 1-. 



\l + a 2 e e 



1 ~ / 

l s ~ l\l + a 2 e d \ , 
1ft ~ + a 2 e e \ , 





d 2 








k~l 


d 2 








l T ~ I 


d 2 







(19) 
(20) 
(21) 



a a€d a n b 2 f K 1 + a2e d)\ , no , 

912 "-rr^' * 23 ^ ed ' 913 ^ P ^^-— d — j. (22) 

where % = (c + abed), and since c u <§C c d e? the quark mixing angles are well approximated 
by the "down" sector. We will not repeat the whole analysis that can be found in [S] but just 
quote the main steps that demonstrate how the horizontal and vertical hierarchy cooperate. 
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- First from lf> ~ Z T , one should have |e e | ~ |e^| as already discussed. 

- Second, one notes that, ignoring the corrections, the second generation Yukawa are unified 

- However, the RG invariant relation s\2 — \Jrndfm s tells us that |a 2 e^| ~ 1, so that a 
should actually be quite large, \a\ ~ |ed| -1 ^ 2 (— 7, see later). 

This large value of a generates the right Cabibbo angle, but also enters the correction of 
the eigenvalues in the first two generations of D, E sectors. 

- Indeed right splitting of s,^ Yukawa is produced by the factors Aj^ = |1 + a 2 ed, e \ once 
the sign (the phases) of a 2 e e and a 2 ed are nearly opposite, a 2 e e ~ —a 2 ed ~ e 0,5 \ so that 
\A e \ ~ 2, \Ad\ ^ 0.5. As discussed above, for the first generation the split goes (correctly) 
in the opposite direction, realizing the small flavour seesaw. 

So, in this ansatz, the largeness of a, following from the Cabibbo angle, explains also the 
deviation from exactly hierarchical family masses of e, d, fi, s. This was already observed 
ini. 

- From other known mass ratios one can then determine \e ej d\ ~ 1/55, a ~ 7.5, |e u | ~ 1/600, 
d 2 ~ 0.5, tan/? ~ 7.3. Then further insight can be gained by checking the CP phase and 
the other mixing angles: from these it turns out that also b is quite large, like a, i.e. 
\b\ ~ 5. Also, one finds two branches: in one the complex b and c are almost orthogonal, 
in the other they are opposite. Finally, there is a flat direction, a combination of | e e / | 
with the common phase of 6, c. 

Summing up, the ansatz is solved successfully and shows a one parameter family of solutions, 
with two branches. The precise leading-order solution for all parameters is not conclusive 
since next-to-leading corrections and contributions to CKM from the "up" sector are 10% 
in magnitude, larger than experimental errors; hence a numeric fit is necessary and will be 
described below. Nevertheless the ansatz is capable of accounting for all the charged fermion 
masses and mixings, unifying the determinants and giving an additional mechanism that links 
the deviations A e< d with the Cabibbo angle 9\2- 

Before describing a model that realizes this ansatz, we would like to comment on the solution 
found from the charged fermions so far. Collecting the three p n vectors: 

pi = (1, 7, 5) , P2 ~ (0, 1, 0.2-1.5) , p 3 ~ (0, 0, 0.7) , (23) 

we observe that they tend to lie in the 2-3 plane. Exact planarity is broken to order 1/8 ~ 
as can be seen from the first 1 with respect to the modulus of its vector \p\ \ ~ 8.5. The reason 
for this is to be tracked in the magnitude of the quark mixing angles that require a and b to be 
large while c and d are of order one. 5 

We find this pattern of quasi-planar triplets a nice hint for the realization of a flavour sector 
of the theory, where the p n vectors are generated by the VEV of three scalar fields that are 
triplets of the flavour SU(3). The potential and the radiative corrections to it lead to their 
mutual disorientation, along the liens of |121 I13j . Incidentally, we also note that while we 
started from complex triplets, fitting fermion masses and mixings indicates that the only two 
complex parameters b, c have a common phase. This means that the three triplets VEV can 
be taken completely real. The CP violation is then produced just by the complex ej and the 

5 We observe also that if one may be concerned with the moduli of p n not being equal and of order 1, (they 
are in fact \pi \ ~ 8, \p2 \ 2^ 1.5, \pg\ ~ 0.7), they can be brought to be in the range 1-2 by rescaling all the e/ by a 
factor of 10. With this normalization the hierarchy of eigenvalues is due partially to the e/ and partially to the 
p n vectors having hierarchic entries, while the normalization that we adopted above corresponds to an eigenvalue 
hierarchy that comes only from the e/'s. The choice may be varied depending on the theoretical realization. 
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coupling constants. This is most welcome, since in usual VEV configurations, the triplets tend 
to repel or attract each other and it is hard to obtain complex relative angles. We leave this 
analysis and the realization of the flavour sector for future work. 

In the next section we complete the description of the model in the SO(10) side, where 
the rank-one projectors are mixed with the SO(10) breaking VEVs that generate the hierarchy 
parameters. 

4 SO (10) universal seesaw and flavour 

Since the different e u ^,e hierarchies correspond to a broken SO(10), it is natural to generate 
them at SO(10) breaking, i.e. at the GUT scale, and then transfer them in some way to the 
light fermions-higgs Yukawa couplings. At the GUT scale we should also mix them with the 
flavour triplets VEVs p n mentioned above. We describe here a way to implement this program, 
along the lines of |S1 IB] • 

4.1 Small parameters from SO (10) breaking 

First we recall that among the different schemes of SO(10) breaking, a very successful and 
economical one is that implementing the so called missing- VEV mechanism [Uj, that while 
breaking to the SM gauge group, automatically stabilizes the hierarchy of the light doublets 
with respect to the GUT scale. Among other fields, in this mechanism there are three multiplets 
in the 45 representation: 45x, 45#, 45#l- Altogether they break SO(10) to the SM group, 
therefore their VEVs are, in Pati-Salam notation, only allowed along the two directions (15, 1, 1) 
and (1, 1, 3). It turns out that (45 R ) G (1, 1, 3), (45 BL ) G (15, 1, 1) and (45*) G (15, 1, 1) + 
(1,1,3). 

All these VEVs can be used in a ratio to some higher scale (e.g. M ~ 10 17 GeV) to introduce 
small parameters into the game. However the last one, the VEV of 45*, is of particular interest 
since it depends generically on two coefficients, that define its mixture of (15, 1, 1) and (1, 1, 3). 
When the small ratio (45x)/M is projected onto the fermion multiplets and it decomposes in 
four small parameters e u €d e e e UD , we have that only two of these are independent, and they 
turn out [5| to be constrained by the following relations: 

e e = -e<i - 2e u , e VD = -2e d - e u . (24) 

The first equation tells us that once \e u \ <C \td\-, the equality in the D,E sectors |e e | ~ |e^| 
is predicted, nicely reproducing the pattern observed in nature. Moreover, the same relation 
implies also bottom-tau unification, with a deviation of e u /e e ~ 10%, that is indeed observed. 

The second equation gives a prediction on the Yukawa hierarchy in the neutrino sector, and 
implies that it is milder than that of the leptons, e VD ~ 1/25. 

Therefore without introducing new fields in the theory, we found that the 45x generates 
the right pattern of small hierarchy parameters e u , e^, e e , e UD . For this reason we will assume 
below that the flavour hierarchies are generated as effective operators by powers of (45x)/M. 

4.2 SO (10) universal seesaw 

To generate the Yukawa couplings as the correct effective operators, it is necessary to introduce 
new (super)heavy fermion multiplets, that act as messengers. Let us start with three families 
of heavy vectorlike fermions 16', 16j. To generate the light Yukawa couplings they should 
couple both to the light fermions and to the Higgs field, that we take to be the usual 10 = 
T(6, 1, + 2, 2). This contains the MSSM higgses in <j). The three families of light fermions 
then receive an effective Yukawa coupling matrix from the mass matrix of heavy vectorlike 
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Figure 3: Universal seesaw mechanism using 45 r. 



fermions 16^, 16^, via a "universal" seesaw mechanism 9 . In detail, the superpotential contains 

W = 7 10 16*16^ + i 4h R 16' 16^ + 16, M^l^ (25) 

where 10 is the fundamental Higgs and is an adjoint one whose role will be clear in a 
moment. Notice that by SU(3) flavour symmetry a direct Yukawa coupling 161610 is forbidden. 
The flavour structure is encoded in Ad heavy, while the first two terms are flavour universal. Near 
at the GUT scale, the higgses develop a VEV, so that the heavy fermions 16' decouple and the 
light fermions 16 acquire an effective Yukawa coupling with 10 that is approximately given by 
the "seesaw" formula: 

l = 17 '(45 R )M- e 1 avy , (26) 

Therefore the Yukawa couplings of the light fermions I are proportional to the inverse of the 
heavy fermion masses, that in turn should exhibit an inverted hierarchy pattern, with the 
lightest particle being the heavy correspondent of the "top" . 

To make explicit the universal seesaw following from (|25ft . we decompose in Pati-Salam 
notation the fermion fields as: 16* = f(4, 2, 1) + f ic (4, 1, 2) and 16^ = ^(4, 2, 1) + Ff (4, 1, 2), 
16j = ^(4,2,1) + Fi(4, 1,2). Note that / = q,l and T = q',l' are isospin doublets and 
f c = u c , d c , e c , u c , F = u',d',e',i>', are isospin singlets. With this decomposition the couplings 
can be illustrated as follows: 

f c F c T c 
f ( 7 (/> \ 
F M R M f 
T \ 70 M T J 

where Mr is the scale of {45r) and Mh eavy is decomposed in Mp and Mjr , respectively the 
mass matrices of heavy isosinglets and isodoublets (unmixed) . 

Only the F, F c participate in the decoupling process and in each sector the Yukawa matrix 
is proportional to the inverse of Mp alone: 

i } = 1 M R M- x , f = u,d,e,u D . (27) 

It is important to note that the j-T c entry is zero because {45r} is in the "right" direction. 
This shows three remarkable features of this universal seesaw using 45r G (1, 1, 3): 1) only the 
isosinglets F c participate in the seesaw mixing with / c , so that their mass matrix is directly 
reproduced in the light Yukawa; this avoids a mixing with Mjr that could spoil the exact 
hierarchies; 2) the LLLL (dominant) part of the D=5 proton decay is automatically eliminated, 
because (1,1,3) does not couple the light and heavy isodoublets [H|ll0j: 3) finally, again because 
the heavy isodoublets do not enter the seesaw, the squark and sleptons mass matrices are 
automatically aligned with the (square of the) Yukawa matrices, and this eliminates problems 
of flavour changing neutral currents. We conclude that here the universal seesaw via the SU(2)^ 
breaking automatically allows the correct mass generation, suppresses proton decay and SUSY 
flavour problems. 
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4.3 Entangling flavour and SO(10) higgses 

We now come to the generation of the flavour structure in Mheavy . We have already introduced 
the three flavour projectors P n that should be coupled to different powers of the 45x field. 
Since however from l|2fi|). (j27|) the heavy mass matrix is the inverse of the Yukawa couplings, we 
must take advantage of a property of combinations of rank-one projectors: the inverse of such 
a combination is a combination of new projectors, with inverse coefficients: 

(aiA + a 2 P 2 + asPs)' 1 = ( — Qi + — Q2 + —Qs) , (28) 

where the Q n are the "reciprocal" of the P n projectors: if P n are parametrized with flavour 
triplets as P n = p n Pn, then Q n = (/ n </* are given by three new triplets q n that are in fact the 
reciprocal of the p n ones: q n = \e nmT p m Ap r /D p , with D p =p 1 - (p 2 Ap 3 ). 

Therefore one should generate a Mheavy as a combination of three rank-one projectors Q n 
multiplied by powers of 45x/M. To have hierarchic eigenvalues these should be successive 
powers, i.e.: 



M M 2 

We do not specify at this point the effective realization of this operator, that will be addressed 
later on. We only observe that such effective operators may involve integer powers (positive or 
negative) of fields, and therefore one may have in front a generic power of it: 



M heavy oc {^fy n [Q t + ^Q 2 + Q 3 ] , (29) 



It is clear that when at SO(10) breaking (45x)/M generates the small parameters e u , e^, e e , 
e UD , one will receive effective Yukawa couplings in the form: 



n-2 



If = ^—[e}P 1 + e f P 2 + P 3 ] f = u,d,e, u D (30) 



e 

I 

The interesting cases are n = 0,1,2, that correspond to the three cases of figure ^ 



a) with n = 0, the coefficients are l,e^ ,e^ 2 , and the GUT Yukawa couplings of the first 
generation are approximately unified. This requires very low tan/3 (< 2); 

b) with n = 1, the coefficients are ef, 1, ej 1 , we have an approximate unification of the second 
generation, at intermediate tan/3 (~ 10); 

c) with n = 2, the coefficients are e^e'A, and large tan/3 (> 50), one reaches unification of 
the third generation, as in the known case of t-b-T unification. 

The second case is the most interesting since there is optimal agreement with the data, and 
as discussed in the present work, the determinants are unified. For the first case see jS], for a 
unified treatment see 

Sticking to the interesting case of n=l (b), we end this section by providing the complete 
realization with renormalizable operators. Such a realization is not unique, so we display two 
different ways to achieve the same result. 

Every realization of nonrenormalizable operators like ()29|) requires one to introduce addi- 
tional fields. For example one can introduce four more vector- like fermion multiplets that are 
flavour singlets: 16", I62, I63 and 16'" (and their conjugates); as an other example one can use 
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in place of three triplets q n , two triplets gi i2 plus one sextet X3> and then introduce three vector- 
like fermion fields: two flavour singlets 16" 2 and one triplet 16'"* (and their SO(10)xSU(3) 
conjugates). The coupling are arranged as follows in these two cases: 





16* 


16' 


16',' 
u 1 


16^' 


16',' 


16"' 




16* 


16' 


16'/ 




16"' 2 


-LO 




1 n 










1 fii 




1 n 








16' 


45 D 




a j 

HI 


( i2 




n 


16' 








'h 


45 v 








45 v 








16? 




?i 


45 X 






16 2 ' 








M 






16 2 ' 




55 




M 




163 




4 








M 












M 


16"' 











M 


45 X 





In the first example q n define exactly the three reciprocal triplets mentioned above, while in 
the second the Qi,2,3 projectors will be given by Q\ oc q%q\, Q2 oc q 2 q\ and Q3 oc X3- 

One final remark is in order: in these examples some of the couplings are forbidden. We do 
not address here the way to motivate this, that can be done by means of additional symmetries, 
discrete or continuous. Indeed the way that heavy higgses and flavour triplets/sextets are 
coupled in these examples (and e.g. in (J2HJ)) points towards a further symmetry that explains 
the charges or couplings much in the Froggatt-Nielsen spirit. 



5 Neutrinos 

In this section we analyze how the neutrino sector can be accommodated in the present frame- 
work of flavour triplets coupled with SO (10) heavy Higgs fields. From the mechanism described 
above, the Dirac neutrino masses are unified together with the other fermions and thus result in 
the range 10-1000MeV. Therefore we need to introduce also a Majorana mass for right-handed 
neutrinos, and generate the light neutrino mass by canonical seesaw [llj . Then given that in 
the model the flavour structure is encoded in the three P n projectors, we assume that also the 
Majorana masses are built with the same projectors. 

The Majorana mass for the RH neutrinos that belong to the light 16 multiplet is generated 
via an universal seesaw similar to that of Yukawa couplings, this time via a singlet messenger. 
We do not dwell in the details that may be found in |Hj, and we quote the result, that also 
the RH neutrino mass M UR will be built with the same flavour projectors P n and with a new 
hierarchic parameter e UR . The coefficients however are generic, since the three new coupling 
constants can no more be absorbed in the P n . We choose to make explicit the e VR and add a 
further complex parameter a. 



l VR oc 



Pi + evn p 2 + —Ps 
a 



(31) 



where for simplicity here and in the following we prefer to use the dimensionless matrices 
Ivr,vl = M VR ,v L /v sin/3, with the electroweak v = 174 GeV. 

A mass matrix for the left-handed neutrinos is then generated by canonical seesaw as l v = 
lv D lv R lu D , between the neutrino yukawa (|17|) and the RH mass matrix (|31j) . The result can 
be found by exploiting an other interesting property of rank-one combinations: indeed, the 
following "seesaw" relation holds: 

(a x A + a 2 A + a 3 P 3 ) (fa P\ + ft P 2 + ft A) _1 (t*i A + a 2 A + «3 A) = 

= U Pi+ s ft+ s4 (32) 
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i.e. the seesaw acts on the coefficients only, but does not change the flavour projectors P n . 
Hence, we get that l VL is built again as a combination of the same flavour projectors: 



(33) 



Using the explicit form (|17|). we can rewrite: 



oc 



\e u /A u 
\ae u /A l 
\be v /A v 



\ae u /A u 
1 
P 



\be v /A v \ 
P 

P 2 + (T 



with A v = (1 + a 2 e u ), p = (c + abe v )/A v , a = (d 2 a + (b - ac) 2 e 2 ,) / A u . 

In this matrix, we only have two free complex parameters e u and a, 6 therefore the model 
is quite constrained since it must explain the neutrino mass hierarchy and the three mixing 
angles. 

Unlike the case of charged fermions, an exact analytical solution of this matrix is not possible, 
due to the large neutrino mixings and to the small but important uncertainties in the parameters 
determined from the charged fermions. Nevertheless it is possible to see that neutrinos are 
predicted to be hierarchic, non degenerate, and to derive a prediction for the #13 mixing angle. 
The analysis goes as follows: 

- Non-degenerate neutrinos with m U 2/m u3 ~ 1/5 follow from the fact that a is large, so the 
ratio of the 1-1 with the 1-2, 1-3 entries is small ~ 1/7. 

- Non-zero 9 e3 follows from a > b, so that the 1-2 entry is larger than the 1-3 one. An 
analytical estimate of the angles from the 2,3 sector can be given first in the form of 
a relation between them, using the parameters found from the charged sector and after 
imposing maximal Q a tm- 



<?e3 



m u2 
m u3 



1 + Re 



16 



1 - Re 1 



with R 



b 

a 



qquark 
'13 



1 



rfluark 
'12 



1.6 ~ 2.1 



(34) 



where 5 = arg p = arg [(c+abe u )/(l + a 2 e u )}. After imposing the correct neutrino hierarchy 
and solving the model for e v , we can estimate also: 



'12 



30° 



16 



; 13 



6 C 



l + Re v 



(35) 



Hence, choosing e purely imaginary or just slightly negative, we see that the right 9 so i ~ 
32° is accomodated, with a prediction for 9 e3 ~ 4°. 

The precise numbers are sensible to the quark mixing angles and CP-phase, as well as 
to the neutrino hierarchy. In practice 9 e3 can be lowered to reach also 9 e3 ~ 0.5° by 
stretching within la the angles and the neutrino hierarchy. 

- Finally, a new flat direction emerges, along which the angles and hierarchy do not vary. 
It directly corresponds to the leptonic CP violation phase, that may be predicted once 
angles and hierarchy will be known to better accuracy. Also a majorana phase varies 
along this flat direction. 

Summing up, the framework naturally predicts hierarchic neutrinos with nonzero 9 e3 in the 
la range, and there is a useful conspiracy from the charged sector to allow the right neutrino 
pattern. A new flat direction emerges as a combination of the four real parameters, to be added 
to the one found in the charged leptons sector. Therefore the model effectively takes advantage 
of just 13 of the 15 real parameters to reproduce the known data. 

6 The common phase of b, c, that was a flat direction of the charged fermions can be reabsorbed in the phase 
of a. 
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Table 1: Best fit with M SU sy 
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Table 2: Best fit with M SU sy = 1 TeV total x 2 = 0.14 per d.o.f. 



5.1 Numeric Fit and Results 

As mentioned above, we performed in |Sj a numerical analysis, taking into account the exact 
renormalization factors. We performed a best fit of neutrinos together with charged fermions. 
The # e 3 neutrino mixing, that we have shown to be predicted inside its allowed la range, is not 
used as a fit condition, but rather displayed as a model prediction. The fit was performed using 
the 15 parameters (11 charged fermions, 4 neutrino) against all known data (20 tests). 

The best fit results are shown in tables ^ and |5J the overall x 2 = 0-35 per d.o.f. is quite 
good, and the model can account for all experimental constraints, albeit having a bit low m u 
and a bit large m s (la level) as implied by the predictions described in section |21 

The fit results confirm the analytical study: 

- Apart from m u , m s that deviate at most 1-a, all data are fitted without tension (< 0.3o"). 

- The flat direction and the two branches correctly emerge from the fit. In the second 
branch all complex phases are almost aligned, hinting for a model with complex aligned 
ej's and real triplets, with a reduction of the total number of real parameters to 10. 

- Neutrinos are predicted as hierarchic and non degenerate, and the mixing from terrestrial 
neutrinos is nonzero: 6 e ^ ~ 2.5°, 4.9°. 
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- Neutrino come out hierarchic and non degenerate m v ~ (0.0045, 0.01, 0.048) eV, and 
right-handed neutrino masses are of the order (10 7 , 10 9 , 10 11 ) GeV. 

6 Discussion 

In this work we discussed the structure of the Yukawa matrices as renormalized at a common 
high scale, and in the context of MSSM GUT, noticed the unification of their determinants. 
This new unification translates in precise predictions in terms of charged fermion masses: in 
particular it gives a prediction for mdm s in la accordance with the experimental ranges, and 
a prediction of tan j3 ~ 7 that lies in the currently favoured range. In addition this mechanism 
introduces a constraint that may explain why the first family fermions are light when the third 
ones are heavy, in a kind of flavour seesaw. 

In fact, as opposed to traditional unification schemes, this idea does not unify some combi- 
nation of yukawa matrices at GUT scale, but unifies approximately the yukawa eigenvalues of 
the second family. 

Examining the mass matrices in more detail we presented an ansatz, similar to the ones 
put forward in in [3J , that realizes the unification of determinants by building the yukawa 
matrices as hierarchical combinations of three generic rank-one projectors in flavour space. This 
construction is successful in explaining the deviations from exact mass hierarchy in the D, E 
sectors, linking them to the Cabibbo mixing. 

We reviewed how these ideas can be realized in a predictive model of SU(3) flavour symmetry 
in the context of SO(10) SUSY GUT, along the lines of jS]. In this model the mass matrices 
can be analyzed with a leading-order analytical approach and with a complete numerical fit. 
All fermion masses and mixings are accomodated at la level, including neutrinos, for which the 
analysis and the numerical fit suggest direct and hierarchical neutrinos with non zero 9 e ^ = 1.7°- 
4.3° or 0.5°-7°, in two branches. 

The model uses the tools that are already available in the context of SO(10) SUSY GUT: 
a "universal seesaw" mechanism and the higgses present in the missing VEV mechanism. The 
universal seesaw is used to to transfer the mass matrices from a set of superheavy fermions 
living near GUT to the ordinary light ones, the higgses in the 45 representation that realize the 
missing VEV mechanism automatically generate the correct horizontal hierarchies e u <C e e ,d- 
Then the mass matrices are generated by coupling one of these higgses with three flavour triplets 
scalar fields, that define three rank-one projectors. This model realization via universal seesaw 
efficiently suppresses proton decay and SUSY FCNC, and maintains almost pure MSSM content 
below the GUT scale. 

In the flavour sector, the solution emerging from the data-fitting points to quasi-aligned or 
quasi-planar triplets, that seem to suggest a precise mechanism of SU(3) breaking along VEVs 
without introducing ad-hoc breakings of flavour, as described for example in |12| I13j. 

We conclude that the unification of Yukawa determinants at the GUT scale nicely fits into 
the framework of SO(10) Grand Unification with horizontal symmetries, where it gives a hint 
for building the flavour sector of the theory, while still giving model independent predictions 
for the fermion masses and the the tan (3 MSSM parameter. 
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A Appendix: MSSM renormalization 



In the MSSM, fermion masses and mixing angles can be denned as follows in terms of quantities 
at the GUT scale: 

m u = l u v sin (3 rj uds R u Bf m d = l d v cos (3 r] uds R d m e = l e v cos (3 r] efl R e 

m c = l c v sin (3 r] c R u B? m s = l s v cos (3 r] uds R d = l^v cos (3 r] efl R e 

rat = hv sin (3 7] t R u Bf mi, = l^v cos /3 r]bR d B t m T = l T v cos (3rj T R e (36) 

Gl2\z = G\2i #13 U = #12 B t 1 , #23 1 Z = #12 B t 1 , (37) 

where the factors R u ,d,e account for the running induced by the MSSM gauge sector, from Mqut 
to nit; the factor B t accounts for the running induced by the large If, and n eiiudsbc complete 
the QCD+QED running from m t down to 2GeV for u, d, s or to the respective masses for b, c. 
From [3], assuming a^\z = 0.118 and M$usy = rnu we fi n d 

7^ = 1.74, 7/c = 2.11, 776 = 1.52, ^ = 1.04, ^ = 1.06, n t = 1 . (38) 

Then we calculate 



R u = 3.45 , R d = 3.36 , R e = 1.51 , R VD = 1.39 , (39) 

Finally, the factor Bt is a function of It- We obtain -E>t|; t= o.5 = 0.901, i?t |zt=3.5 = 0.670. Near 
lower /t, as in the case of moderate or large tan/3, the dependence on It is stiffer, and in power 
6 it allows a quite precise dermination of It from the experimental valoue of mt- 

The dependence on a^\z of all these coefficients affects mainly the lighest quarks via rj uds , 
7] c (10%), all the others have uncertainty of the order of 1-3%. There is a stronger dependence 
on the supersymmetry breaking scale. 
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